Universal infrared conductivity of graphite 
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The conductivity of graphite is analyticaUy evaluated in the range of 0.1-1.5 eV, where the electron 
relaxation processes can be neglected, and the low energy excitations at the "Dirac" points are most 
essential. The value of conductivity calculated per one graphite layer is close to the universal 
conductivity of graphene. The features of the conductivity are explained in terms of singularities of 
the electron dispersion in graphite. 
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Since the pioneering experimental investigations of a 
single atomic layer of graphite (graphene)^'^ , its prop- 
erties attract much attention. Among them, the optical 
response is of particular interest. Recently the transmit- 
tance of light throw the graphene monolayer has been 
measured^Ti^. The transmittance 

T = 1 — TTQ! 

was found to be frequency independent in a broad range 
of photon energy. The result of the experiments is re- 
markable because it involves the fine structure constant 
a. It was discovered that the real part of the optical 
conductance of graphene takes the universal value 
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which does not depend on any parameters of graphene. 
This value agrees perfectly with the calculations^i^ ig- 
noring the Coulomb interactions between electrons. The 
agreement shows that the poorly screened Coulomb in- 
teraction does not play any role in graphene for infrared 
photon frequencies^'^. 

The intermediate place between 2d graphene and 3d 
semiconductors belongs to multilayer grapheneaiS, and 
graphite, which have a layered structure with the inter- 
layer distance Cq = 3.35A much larger than the nearest- 
neighbor distance ao = 1A2A in the layer. In the study of 
graphiteii, it was found that its optical conductivity per 
one layer is very closed to the universal conductivity of 
graphene and has evident peculiarities. The analytic cal- 
culation of the in-plane optical response of graphite done 
previouslyiS has ignored coupling between layers and no 
peculiarities have appeared for the infrared region. 

In the present paper, we evaluate analyticaly the con- 
ductance of graphite in the infrared region of the photon 
frequencies. It is known that the low energy electron ex- 
citations in graphene can be described very well with 
the Slonczewski-Weiss-McClure theoryi^. The largest 
parameter of the theory, 70 = 3.1 eV— , describes the 
electron dispersion for in-layer directions k. If the photon 
energy is less than 70 , we can use the linear expansion of 
the in-layer hopping term in the Hamiltonian and intro- 
duce the constant velocity parameter v = 10^ cm/s. The 



second parameter in the rang is the interlayer hopping 
71 of the order of 0.4 eV which is known from experi- 
ments on bilayer graphen o^^i^^ . The parameters 73 and 
74 give the corrections of the order of 10% to the in-layer 
velocity v. The electron-hole overlap of the order of 0.02 
eV is determined by parameters 72 and 75 (see Fig. [5]). 
Therefore, for the photon frequencies larger than 0.1 eV, 
we can neglect the terms with 72 and 75. Calculating 
such the integral property as conductivity in the region 
of the infrared frequencies between 0.1 eV and 1.5 eV, 
we can, first, neglect the small parameters of the theory 
and, second, use the fc-expansion of the in-layer hopping 
term. Our results have the evident analytic form. 

In this approximation, the effective Hamiltonian writes 
near the K-G-H lines of the Brillouin zone in the simple 
form 



«+ liz) 



iJ(k) = 





fc_ 

j{z) fc_ 

fc+ 



(1) 



determined only by two constants. One is v — lO^cm/s 
included in the definition of the in-plane momentum com- 
ponents, k± = v{^ikx — ky), and another is the inter-layer 
interaction 71 involved in the function 7(2;) = 271 cos 2. 
The momentum component z = kzCo is limited by the 
Brillouin half-zone, < z < 7r/2 in relative units. 
The corresponding eigenenergies are 
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On the K-G-H lines, k — 0, these equations determine 
two bands £1^4 = ^^iz) and two degenerate (electron 
and hole) bands with the energy £2,3 = 0. We have to 
emphasize that this degeneracy results from C31, symme- 
try on the K-G-H line. 

In order to calculate the conductivity, we use the gen- 
eral expression 
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FIG. 1: The dispersion of the low energy electron bands in 
graphite. 



has the matrix elements 

V23 = 2i{ez - e2){-kxe. 

vi2 = 2(£i + e2){Kex + kyey)]/NiN2 
vi4 = 2i{ei - ei){-Key 



kyex)]/N2N3 



kyex)]/NiNi 



where are the unit vectors directed along the coor- 
dinate axes. For the real part of conductivity, the in- 
tegration in Eq. ^ is easily taken at zero temper- 
atures T = in cylindrical coordinates {kz,k,(p) over 
the angle and over k with the help of the i5-function, 
(lo — X + iv)^^ — > — Z7r(5(a; — x). One obtains for con- 
tributions of the transitions between the corresponding 
valence and conduction bands into the diagonal compo- 
nents of conductivity (off-diagonal ones equal zero) the 
following integrals over z = kz/c^. 
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valid in the collisionless limit 3> i^, where u is the 
collision rate of the carriers, /(e) — [exp(^^^) — 1]^-'^ is 
the Fermi-Dirac distribution function, and the integral is 
over the Brillouin zone. 

Here, the first term is the Drude-Boltzmann conductiv- 
ity negligible for frequencies larger than the electron-hole 
overlap. The second term represents the optical inter- 
band transitions of electrons from the valence 2,4 to con- 
ductive 1,3 bands. The real part of the interband contri- 
butions into conductivity arises from the bypass around 
the pole at en(k) — £m(k) = ±w. The imaginary part is 
given by the principal value of the integral. 

The velocity operator 

_ 9g(k) 

^ ^ ak 

near the K-G-H lines is determined by the Hamiltonian 
((!]) . The corresponding matrix elements should be calcu- 
lated in the representation, where the Hamiltonian has a 
diagonal form. The operator transforming the Hamilto- 
nian to this form can be written as follows 
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where 7(2;) = 271 cosz and 0{x) is the step function. 

It is evident from Eqs. ^ (see also Fig. [2]) that the 
conductivity (T23 tends to e'^ /AHcq at the low frequencies 
Lu <^ 271, whereas other contributions go to zero in the 
limit of low frequencies. At larger frequencies uj >• 271, 
the total conductivity (the sum of (T23 and 1741) tends 
again to e'^/Ahco- Therefore, ctq — e^/4?ico can be con- 
sidered as the universal conductivity of graphite, where 

/Ah is the conductivity of monolayer graphene and the 
factor 1/co is the number of the layers per the length unit 
in the z-direction of graphite. Integrating in Eqs. (jSj, we 
get finally 
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where = 2{el 
locity operator 



fc^) . In this representation, the ve- 
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FIG. 2: The real part of the graphite conductivity per layer 
(in units of /Ah) versus the frequency (in units of 271 — 
0.84 eV); the experimental dataii are shown in the solid line, 
results of the present theory in the dashed line. The insert 
shows the contributions of various electron transitions. 



and shown in Fig. |3l Here, the peculiarity looks like 
a threshold at = 271 and it is more clearly marked 
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FIG. 3: The imaginary part of the graphite conductivity per 
layer (in units of /Ah) versus the frequency (in units of 
271 = 0.84 eV). 
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where t = w/271 and zi = arccos(i/2). 

The peculiarity as a kink can be seen in Fig. [2] The 
expression ([5]) shows that this kink is located at cj = 
271. Taking into account the kink position ut — 0.84 eV 
determined experimentally, the value of 71 = 0.42 eV is 
found in excellent agreement with experiments on bilayer 
graphene. 

The contributions of the electron interband transitions 
into the imaginary part of conductivity can be integrated 
over k at the zero temperature. The results are obtained 
in the form of integrals over kz 
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in comparison with the kink in the real conductivity. 
Both peculiarities result due to the electron transitions 
between the bands 2 — >■ 1 and 4 — )■ 3. We should em- 
phasize that the peculiarities become broader with the 
temperatures and the collision processes included. 

So far the in-layer conductivity was considered. The 
estimate of the inter-layer conductivity can be also done. 
Since the conductivity is determined by the ratio of the 
corresponding velocities squared, we have to write 

Vz = ~ 71 Co sm(/ezCo) . 

Then, integrating over fc^, we get 

(Jz/ao^ {liCo/hv)^/2^Qm. 

In conclusion, our calculations reveals that the optical 
conductance of graphite can be estimated for frequencies 
between 0.1 and 1.5 eV multiplying the graphene con- 
ductivity /Ah by the number of the layers 1/co per the 
length unit. The Drude-Boltzmann contribution is essen- 
tial at lower frequencies, whereas others interband transi- 
tions, e.g. at the M point of the Brillouin zone contribute 
into the conductivity at higher frequencies. The similar 
estimate are applicable for other graphite materials such 
as nanoribbons. The kink in the real part of conductiv- 
ity and the threshold in the imaginary part appear at 
the frequency uj — 271 determined by the interlayer cou- 
pling. The sharpness of the features are smeared with 
the relaxation processes and temperatures included. 
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